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ABSTRACT 

We investigate the dynamical status of the low-mass globular cluster Palomar 13 by 
means of iV-body computations to test whether its unusually high mass-to-light ratio 
of about 40 and its peculiarly shallow surface density profile can be caused by tidal 
shocking. Alternatively, we test - by varying the assumed proper motion - if the 
orbital phase of Palomar 13 within its orbit about the Milky Way can influence its 
ap pearance and thus ma y be the origin of these peculiarities, as has been suggested 
by iKupper et al.l (|2010bf) . We find that, of these two scenarios, only the latter can 
explain the observed mass-to-light ratio and surface density profile. We note, however, 
that the particular orbit that best reproduces those observed parameters has a proper 
motion inconsistent with the available literature value. We discuss this discrepancy and 
suggest that it may be caused by an underestimation of the observational uncertainties 
in the proper motion determination. We demonstrate that Palomar 13 is most likely 
near apogalacticon, which makes the cluster appear supervirial and blown-up due to 
orbital compression of its tidal debris. Since the satellites of the Milky Way are on 
average closer to apo- than perigalacticon, their internal dynamics may be influenced 
by the same effect, and we advocate that this needs to be taken into account when 
interpreting their kinematical data. Moreover, we briefly discuss the influence of a 
possible binary population on such measurements. 

Key words: galaxies: kinematics and dynamics - galaxies: star clusters - globular- 
clusters: individual: Palomar 13 - methods: ./V-body simulations 



1 INTRODUCTION 

There are many objects on the sky, especially in the halo of 
the Milky Way (MW) , whose nature is not clear to us. Some 
of those objects are hard to address observationally, and for 
others there is just no conclusive theoretical explanation. 

In fact, simply by looking at a colour-selected sample 
of stars within a region of the sky it is sometimes not easy 
to determine the true extent of a stellar system, mostly 
since it lacks a clear cut-off in its surface density profile. 
The same holds true for the determination of its velocity 
dispersion through a sub-sample of stars with readily mea- 
sured radial velocities. These uncertainties typically result 
in discussions and speculations about a best-fitt ing density 
profile as well as a system's true tidal radius (e.g. King 1966; 
Elson. Fall fc Freeman] 1 19871 : iMcLaughlin fc van der Marell 



20051 ). and a lso about the true mass-to-light r atios of such 



systems (e.g. lKroupalll997l ; [Mieske et al.ll2008h . 
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Some of these uncertainties arise from peculiar sur- 
face density profiles. That is, even though many objects 
in the Milky Way halo are well limited and show a well 
defined surface density profile with a slope of about R~ 4 
in the region of the tidal radius, some objects obey shal- 
low surface density profiles in the outskirts, having slopes 
of about -1 to -2 , like for example the M W globular clus- 



ters Palomar 5 l|Odenkirchen et al I l2003f>. NGC 54 66, M 
15 , M 53, M 30, and NGC 5053 (IChun et al.1 l201Ch. AM 



4 (ICarraro. Zinn fc Moni Bidinll2007h. Whitin g 1 ifcarrarol 
l2009f ). and NGC 1851 (|Qlszewski et al.l l2009t l. The latter 
furthermore seems to be surrounded by a 500 pc halo of 
stars whose origin is unknown up to now. 

Other uncertainties arise from unusual mass-to-light 
(M/L) ratios of some stellar systems. While most globu- 
lar clusters show mass-to-light ratios of 1-2, Ultra-Compact 
Dwarf galaxies (UCDs) have higher M/L by a factor of 
about two, wherea s dwarf spheroidal galaxies even show 
values of up to 10 3 l|Dabringhausen. Hilker fc Kroupall2008l : 
iGeha et al.1l2009h . These differences are usually ascribed to 
different dark matter contents, catastrophic tidal heating 
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CFH T bvlCote et al.l d2002T) . Fitted to the profile are a King tem - 
plate jKing|ll962^ and a KKBH template (iKiipper et al.ll2010bh . 
Given in the figure are the template values for the fitted tidal ra- 
dius for both templates, R^ mg and R^ KBH respectively, and for 
the extra-tidal slope, r\ (eta). The shallow slope at large radii of 
r\ = 1.9 influences the fit of the King template such that it yields 
a significantly larger value for the tidal radius as does the KKBH 
template. 1 arcmin corresponds to about 7 pc at the assumed 
distance of Pal 13. 



by gravitational shocks, a variation of the IMF, tidally re- 
shaped stellar phase-space distribution functions, contam- 
inations from stellar streams i n the MW halo, or alterna- 
tive gravitational theories (e.g. Kroupa 1997, Gilmorc et al.l 
20071. ISimon fc Gehall2007t iMieske et alj|2008l . lAngusll200Sl. 



Niederstc-Ostholt et al. 



2009). 



The low-mass Galactic globular cluster Palomar 13 is a 
stellar system which shows both, an unclear extent due to a 
shallow surface density profile, and a high veloci ty dispersion 
resul t ing in a mass-to- light ratio of about 40 (Sie gel et al.l 
l200ll ; ICote et al.ll2002T ). Further details on this cluster are 
presented in Sec. [5] In this investigation we demonstrate by 
means of iV-body calculations how these observational re- 
sults can be explained without the need for dark matter, 
tidal heating, binaries or changes in the law of gravity. To 
this end we compute models of Palomar 13 on various or- 
bits about the Galaxy that are consistent with its present- 
day distance and radial velocity with respect to the Sun. 
We show how different such a stellar system can appear in 
different phases of its orbit. Details on the models are de- 
scribed in Sec. The results of these computations and the 
mock observations in which we show how this cluster may 
appear when observed with an 8m-class telescope, are shown 
in Sec. 0] Sec. [5] is a short discussion on the plausibility of 
our findings. Finally, in Sec. [6] we give a short summary and 
conclusions. 



2 PALOMAR 13 

Palomar 13 is an old and metal poor Galactic globular 
cluster . From iso c hrone fits to its colour-magnitude dia- 
gram, ICote et all l|2002l ) find Pal 13 to be about 13-14 
Gyr old, and from spectroscopy they derive a metallicity 
of [Fe/H] = —1.9 ± 0.2. Moreover, it is among the faintest 
objects liste d in the Har ris catalogue of Milky Way globu- 
lar clusters (|Harris!ll99r3 ). With an estimated mass of about 
3000 Mq (assuming a mass-to-light ratio of unity ) it is one of 
the le ast massive glo bular clusters of the Galaxy (|C6te et al.l 
2002). Furthermore, ISiegel et all l|200ll ) argue with proper 
motion measurements that Pal 13 is on an inclined, highly 
eccentric orbit about the Milky Way. Thus, it most proba- 
bly was subject to strong tidal disruption during the last few 
Gyr. Indeed, observations show further peculiarities about 
this specific cluster. First, the cluster shows an unusually 
high velocity dispersion and therefore a very high mass-to- 
light ratio, and second, its surface density profile differs from 
the ones of most other Milky Way globular clusters: 



(i) Corresponding to data by ICote et"all ((2002) , which 
we mainly use to compare with our computations, the 
cluster is located at Galactic longitude of i = 87.° 1 and 
Galactic latitude of b = — 42. °7. Its distance from the 
Sun is Rsun = 24.31}'^ kpc, placing it at a distance of 



Roc 



25.3 4 



kpc from the Galactic centre. Pal 13's 



radial velocity was determined by Cote et al. to be 
V r = (24.1 ±0.5) km/s using spectroscopic data of the High 
Resolution Echelle Spectrometer at the Keck telescope. 
In the same investigation, its internal radial velocity 
dispersion was found to be oy = (2.2 ± 0.4) km/s from 
a sample of 21 stars located within the cluster's inner 
2 arcmin. Within their best estimate of Pal 13's tidal 
radius of 26 arcmin, Cote et al. furthermore measured 
an absolute magnitude of My = —3.8 mag. Assuming 
Pal 13 being in virial equilibrium, this would imply a 
mass-to-light ratio of M/L — 40 jl^*. Cote et al. suggest 
that this unusually high velocity dispersion could be the 
consequence of either a catastrophic heating during a recent 
perigalacticon passage or the presence of a dark matter halo. 

(ii) Fig. [T] shows the surface density profile as obtained 
with the 3 . 6 Can ada-France-Hawaii telescope (CFHT) by 
ICote et all (|2002T ). The profile shows a shallow slope 77 ~ 
—2 out to large radii (10 arcmin correspond to about 70 
pc at the distance of Pal 13), markedly different to the 
much steeper slope rj ~ —4 found in most other sur- 
face density profiles of globular c lusters (compare, e.g., 
iMcLaughlin fc van der Marell 120051 ). This large extent of 
Pal 13 can be interpreted in two ways. First, the stellar 
population at large radii can be part of the cluster such 
that Pal 13 would be a very low concentrated c luster with 
a larg e tidal radius of about 26 arcmin [180 pc] l|C6te et alj 
2002), or, second, Pal 13 can be interpreted as having a 
very pronounced tidal debris and the cluster itself having a 
significantly smalle r tidal radius of about 3 arcmin [20 pc] 
(|Siegel et al.ll200"ll ). To check these oppositional proposals 
for consistency with theoretical expectations, we can make 
a first estimate of Pal 13's true tidal radius, R t , using 
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Figure 2. O rbit of Pal 13 , as de termined using the proper motion 
measured bvlSiegel et al.l ll200lf) and the radial velocity measured 
bv lCote et alj<2002l) . referred to as orbit 1, shown in the meriodal 
plane representation, where r depicts the radial distance from the 
Galactic centre within the Galactic disk and z is the height above 
the Galactic disk. Shown are the last 3 Gyr and 0.5 Gyr into 
the future. Also shown is the resulting orbit after setting the 
transverse velocity, Vt, to zero, such that the cluster has only the 
(more precisely determined) radial velocity measured by Cote et 
al. [orbit 2). The black dotted line depicts the orbit with the 
transverse velocity chosen such that it minimizes Pal 13's orbital 
velocity, V or i,, referred to as orbit 3. 

where G is the gravitational constant, M is Pal 13's present- 
day mass, a nd O its angu lar velocity on its orbit about the 
Milky Way (ISpitze j|l987h . This yields 

7?t~43.7pc. (2) 

when we assume that Pal 13 has a mass of about 3000 Mq 
and is on a circular orbit with an orbital velocity of about 
Vorb = 220 km/s at Rgc = 25.3 kpc. In fact, Pal 13 is more 
likely to be on an eccentric orbit and hence may rather have 
a present-day tidal radius of about 50-100 pc when we take 
into account that its true angular velocity is likely to be 
lower than that of a circular orbit. Thus, the theoretically 
expected range for Pal 13's tidal radius does not agree with 
either of the two observational estimates of Cote et al. and 
Siegel et al. 

2.1 Orbit 

The hypothesis that the above peculiarities ar e caused by 
tidal effects is supported by measurements of ISiegel et al.l 
|200ll) . who, using CCD photometry and 40 years older 
photographic plates, find Pal 13's proper motion to be 
fi a cos8 = (2.30 ± 0.26) mas/yr and [is = (0.27 ± 0.25) 
mas/yr. In Fig. [2] the orbit of Pal 13 within the Milky Way 



is shown in red solid lines for the case of combining the 
proper mo t ion of Siegel et al. with the radial velocity of 
ICote et all l|2002t ). Note that we will refer to this orbit as 
orbit 1 throughout the text. The plotted line corre s ponds 
to the orbit integrated within an lAllen fc Santillanl (1991) 
Milky Way potential for the last 3 Gyr and 0.5 Gyr into 
the future. Distance and orbi t al mo tion of the Sun where 
taken from lDehnen fc Binnevl (|l998h . that is, the Sun is lo- 
cated at a Galactocentric distance of x — 8 kpc, the orbital 
velocity of the local standard of rest (LSR) is 220 km/s in 
y-direction, and the Sun moves with respect to the LSR with 
V x = 10.0 km/s, V y = 5.3 km/s and V z = 7.2 km/s. 

As we can see from Fig. [2] corresponding to these mea- 
surements, Pal 13 has a very elliptical orbit with eccentricity 

Rapo Rperi 82.8 11.1 ^ 

Rapo + Rperi 82.8 + 11.1 

where R apo is its apogalactic distance and Rperi its peri- 
galactic distance, respectively. Therefore, with a Galacto- 
centric d istance of 25 . 3 kpc and the proper motion as mea- 
sured bv lSiegel et al.l (|200ll ). Pal 1 3 is today quite clo se to 
its perigalacticon . Because of this, ISiegel et alj l|200lf ) and 
ICote et all l|2002h suggest that the high velocity dispersion 
of Pal 13 and the shallow slope of its surface density profile 
at large radii may well be due to the last pericentre pas- 
sage which may have heated the cluster violently and may 
have caused a rapid expansion and/or an overspilling over 
the tidal boundaries. 

A comprehensive A-body investigation of low-mass 
globular clusters on eccentric orbits however shows that peri- 
centre passages at such great galactocentric di stances barely 
cause violent mass loss or rapid expansion l|Kiipper et al.l 
l2010al ). Furthermore, a follow-up investigation showed that 
the surface density profiles (a) cannot be used to draw con- 
clusions on the (theoretical) tidal radius of a cluster, and 
(b) only show shallow slopes at large radii when the cluste r 
is close to reaching its apogalacticon (|K(ipper et al.ll2010bh . 
The latter is due to the tidal tails of the cluster which get 
stretched and compressed along the orbit. That is, if a clus- 
ter and its tails move from apogalacticon to perigalacticon 
they get accelerated and stretched, whereas from perigalac- 
ticon to apogalacticon they get decelerated and compressed. 

If the shallow slope in the surface den sity profile of 
Pal 13 is indeed due to the effects described in lKiipper et al.l 
l|2010al lbT). then Pal 13 would need to be currently in a po- 
sition close to apogalacticon. This, however, disagrees with 
the position in its orbit [orbit 1) derived from the proper 
motion measured by Siegel et al. and that puts Pal 13 close 
to perigalacticon (Fig. [2]). In an attempt to resolve this con- 
tradiction, we will in the following define two further test 
orbits which differ in proper motion from the Siegel et al. 
estimate. We will then fully integrate the dynamical evolu- 
tion of Pal 13 for all three orbits to compare the results with 
the observed surface brightness profile, luminosity, and ra- 
dial velocity dispersion (Sec. [4]). The underlying assumption 
for this procedure is that the proper motion measurements 
are the most uncertain observed cluster properties. 

We first force the cluster's transverse velocity (proper 
motion) to zero, while keeping the radial velocity. In Fig. [2] 
we can see that this results in an orbit which is less eccentric 
{Rapo = 49.3 kpc, Rperi = 12.5 kpc, e = 0.60) and in which 
the cluster today is closer to apogalacticon. We will refer to 
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this orbit as orbit 2. In addition, we search for the proper 
motion values which minimize Pal 13's 3D orbital velocity, 
Vorb, since this would yield the orbit in which Pal 13 is clos- 
est to apogalacticon. In fact, by doing this we get a more 
eccentric orbit as the cluster falls deeper into the Galactic 
centre (R apo = 38.5 kpc, R per i = 3.5 kpc, e = 0.83). This or- 
bit will be referred to as orbit 3. It is also depicted in Fig. [2l 
with the corresponding values for the proper motion being 
HcCOsS — 0.72 mas/yr and ^ = —1.2 mas/yr. Note that 
these values differ significantly from the values measured by 
ISiegel et al.1 (|200lh by more than 1 mas/yr in each direction. 

For all three orbits, the cluster is coming from perigalac- 
ticon and approaching apogalacticon. But all three predict 
Pal 13 to be in a different orbital phase, p or b, which we here 
define as 



Porb — 



RGC RGC — Rpar 
\Rgc \ Rapo Rper 



(4) 



and which is constructed to be zero in perigalacticon and 
unity in apogalacticon. This factor p or b gives the fraction 
of the radial distance between R peT i and R apo at which a 
cluster is currently located. In this definition Roc is the time 
derivative of the galactocentric radius, which divided by its 
magnitude adds a minus sign to the orbital phase in case 
the cluster is moving from apogalacticon to perigalacticon. 
For circular orbits p OT b is always zero. 

We introduce the orbital phase, p or b, here in addition 
to the orbital eccentricity, e, since it appears crucial for the 
appearance of the observed effects (|Kiipper et al.ll2~010bl ). A 
systematic study on the dependence of the described effects 
on those two parameters will follow in a future investigation. 
For orbit 1 we get p or b = 0.20, for orbit 2 p or b = 0.35, and 
for orbit 3 p or b = 0.62. 



2.2 Surface density profile 

iKtipper et al.l |2010bl ) introduce a template (KKBH) which 
can be used to reliably measure the slope of a surface density 
profile at large radii. The KKBH template reads as follows: 



MR) = k 



R/Rc 



1 + R/Rc 
1 



V 1 + (R/Rc)* \A + (Rt/Rc 
for radii smaller than fiRt, and 

-ri/64 

1 + 



M R ) = M^Rt 



R 

fj-Rt 



(5) 



(6) 



for R ^ fJ,Rt, where A; is a constant, R c gives a core radius, 
7 the core slope inside R c , and Rt a tidal radius (which 
Kiipper et al. name edge radius to avoid confusion with the 
theoretical tidal radius, since those two correlate only under 
certain circumstances). For radii larger than a fraction fi of 
R t the template changes into a power-law with slope r). The 
exponent 64 in /2 (R) causes the template to change abruptly 

into the power- law slope at fiR t . 

KKBH is based on the template of lKinel (| 19621 ) but is 
modified in two steps: first it allows to have a power-law 
cusp in the core, and second it has an additional extra-tidal 



component in the form of a power-law slope. In this way, 
also more concentrated clusters can be represented for 
which the original King template fails, and furthermore the 
cluster profile can be fitted without being influenced by 
a dominant tidal debris. This effect can be seen in Fig. [T] 
where we applied a regular King template fit to the CFHT 
data by Cote et al. and also a fit of KKBlfl. The King 
template gets significantly influenced by the stellar material 
at large radii and yields a tidal radius of more than 15 
arcmin. The KKBH template assigns a power-law slope of 
r\ = 1.91 ± 0.15 to the tidal debris and yields a tidal radius 
of only 1.9 arcmin. 

By computing iV-body models for all three kinds of or- 
bits given above, and fitting the KKBH template in the same 
way to similarly resolved JV-body data, we will try to repro- 
duce this observed surface density slope r\. Furthermore we 
will measure the velocity dispersion and absolute magnitude 
of the computed clusters in the same way as Cote et al. have 
done, and compare these values with the observational ones. 



3 MODELS 

We computed 45 mo dels of Pal 13 using the collisional N- 
body code NBODY6 (|Aarsetbl2003l ) on the GPU computers 
at AlfA Bonn. We set up 15 different cluster configurations 
using the publicly available tool McLustef0 (Kiipper et al., 
in prep.) . We used a tidally truncated Plummer profile where 
we varied the cluster half-mass radius between 4, 6 and 8 
pc, and the initial mass between 3000, 4000, 5000, 7500 and 
10000 Mq, respectively. The mean mass of the cluster stars 
was in all cases about 0.3 Mq, thus the number of objects 
in the computations were a factor of three times higher. 
Each of these 15 clusters was computed for the last 3 Gyr 
on each of the three different orbits (orbit 1-3) mentioned 
in Sec. [2] (see also Fig. [2} . We focus on the last 3 Gyr of 
evolution since we are only interested in the near by tidal 
debris. Using equation 18 of iKiipper et al.l (|2010al ) for the 
mean drift velocity of stars within the tidal tails, 



v c = ±(4GA/fi) 



1/3 



(7) 



where G is the gravitational constant, M is the cluster mass 
and its angular velocity on its orbit about the Milky Way, 
we can estimate the length of the tidal tails after 3 Gyr of 
evolution if the cluster was on a circular orbit. With the 
same values which we used in eq. [5] we get a drift velocity of 
about 0.77 pc/Myr and thus a minimum length of the tails 
of 2.3 kpc in each direction from the cluster. Note that this 
estimate gets complicated through the fact that Pal 13 is 



1 Kiipper et al. recommend using an additional constant back- 
ground, b, for the KKBH fit to allow for more flexibility at large 
radii, but as it turns out, this is only reasonable with highly re- 
solved (e.g., W-body) data. For less well resolved observational 
data with much fewer data points, it is more reasonable to reduce 
the number of fit parameters to a minimum . We the refore set the 
background, b, mentioned in lKiipper et ah (2010b) to zero, and 
in addition fix the break radius parameter fi to 0.5, as was found 
by Kiipper et al. to be the most plausible value. 

2 www . astro . uni-bonn . de/~akuepper/mcluster/mcluster . html 
or www. astro .uni-bonn. de/~webaiub/german/downloads .php 



The Curious Case of Palomar 13 5 



> 
e 




B-V 



Figure 3. Mock colour-magnitude diagram of one of the com- 
puted clusters at the end of 3 Gyr of -/V-body integration, i.e. 
when the stellar population is 13 Gyr old. The y-axis gives the 
apparent V magnitude as would be observed from the distance of 
the Sun. Random errors which grow exponentially with increasing 
magnitude were applied to both apparent magnitudes, my and 
mg (see text). This particular cluster had an initial mass of 5000 
Mq and a half-mass radius of 8 pc. The dashed box shows the 
region in the colour-magnitude diagram in which we define stars 
to be cluster members. 
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Figur e 4. Colour -magnitude diagram of the Bcsancon model 
dRobin et alj|2003l) covering a 1 deg 2 field around the position 
of Pal 13. Within the region which is occupied by Pal 13 in this 
diagram (dashed box) we count about 1000 stars which will pol- 
lute observations when going down to ray = 25 mag. 



most likely not on a circular orbit. Through the acceleration 
and deceleration on an eccentric orbit the cluster-tail system 
gets periodically stretched and compressed. This estimate 
is therefore only a mean value of the length of the tails. 
Anyway, for investigating the vicinity of Pal 13 this timespan 
seems to be sufficient. 

Since we want to produce a realistic CMD of Pal 13 
with the a ppropriate photometric ob servables, we use the 
SSE code ([Hurley. Pols fc Toutll2000l ) in combination with 
McLustee0 to set up evolved stellar populations of 10 Gyr 
age with a metallicity of [Fe/H] = —1.9. T he populations 
are evolved from a canonical Kroupa IMF (|Kroupall200lh 
ranging from 0.08 Mq to 100 Mq, where compact remnants 
are only kept if their kick velocity which gets assigned to 
them by SSE does not exceed Pal 13's present-day escape 
velocity, v e3C , calculated using 



2GM 

Rh 



(8) 



where M is again the cluster mass and Rh is the cluster's 
half-mass radius, respectively. This treatment is a bit arbi- 
trary since we do not dynamically model the first 10 Gyr 
of the cluster's life and the true retention fraction could be 
both, higher or lower. Another simple treatment would be to 
keep all compact remnants which would have a slight effect 
on the observed mass-to-light ratio as the cluster would then 
have more mass than can be seen in stars. But only about 80 
compact remnants get expelled from a 5000 Mq cluster like 
we model here in the way described above, hence we con- 
sider this to be of secondary importance and concentrate on 
the case of low dark mass in the clusters. 

Those evolved clusters we then feed to NBODY6 to 
evolve them further, chemically and dynamically, up to a 
total age of 13 Gyr. In this way we can concentrate on the 
last few Gyr of dynamical evolution of the cluster, which 
are most important for its present-day structure and nearby 
tidal debris. Hence, we save computational time with this 
te chnique. A similar a pproach has been successfully tried 
bv lHurlev et al.1 (|200ll ) for the open cluster M67. 

The stellar evolution of single stars within NBODY6 is 
also calculated with SSE, a consistent treatment of stellar 
evolution throughout the investigation is therefore guaran- 
teed. From NBODY6 we finally extract the luminosities and 
stellar radii of all stars within the calculations to compute 
their effective temperatures and with this th eir colours an d 
magnitudes in the Johnson-Cousins system (jBesselllll990l ). 
We use the algorithm described in iFloweij (|l996l ) to first 
derive the bolometric correction, BC, and the colour index, 
B — V, and with this the absolute magnitude in the V-band, 
My, and in the B-band, Mb- Together with the distance 
information of each star we can then derive the apparent 
magnitudes, my and m.B, respectively, and can apply a re- 
alistic cut-off at a magnitude limit of, e.g., my = 25 mag as 
would be achieved by an 8m-class telescope in a few minutes 
of integration. 

Since any observation obeys statistical and instrumental 
uncertainties, we furthermore apply a Gaussian-distributed 



3 Note that this version of McLuSTER including SSE is also avail- 
able from the given web address. 
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random error, dm to each apparent magnitude, mv and tub, 
which increases with decreasing brightness as 

dm = \/0-02 2 + 0.07 x io.0 - 4x ( m - 25 ). (9) 

This gives a minimum error of 0.02 mag, and an additional 
uncertainty which is of the order of 0.07 mag at 25 mag 
and decreases with increasing brightness. A final, synthetic 
colour-magnitude diagram of one of the clusters is given in 
Fig. 01 

Objects with low surface densities such as tidal debris 
will be largely affected by background/foreground source 
contamination. We therefore have to estimate the number 
of stars which will pollute our mock observations. For this 
purpose we gener ate an artificial ste llar population using the 
Besancon model (|Robin et al.ll2003l ) for a 1 deg 2 field around 
the position of Pal 13 in the same filter set (Fig(4|. Within 
the CMD region occupied by the cluster model we count 
about 1000 stars which can be mistaken as cluster mem- 
bers, corresponding to 0.28 stars/arcmin 2 . That is, Pal 13 
and its tidal debris will only be visible in those places where 
its surface density exceeds this value. In addition to our es- 
timate and for the sake of being conservative, we will also 
discuss our results using the somewhat higher background 
su rface density of 0. 69 stars/arcmin 2 found observationally 



su riace density or u. t>9 
bv lCote et aTT(|2002h . 



4 RESULTS 

From each computation we take a snapshot after 3 Gyr of 
dynamical evolution as seen from the location of the Sun. At 
this point the stellar population is 13 Gyr old, and should 
resemble the stellar population of Pal 13 within the given 
uncertainties. Stellar maps of a 4 deg 2 region around one 
of the clusters (Mo = 5000, Ro — 8 pc) for each of the 
three orbits are shown in the left panels of Fig. 1 51 71 In these 
figures each dot represents a star above my = 25 mag. A 
background of 0.28 stars/arcmin 2 was added with random 
positions. By comparing the figures we see that the orbit 
with the Siegel et al. proper motion (orbit 1, Fig. [5]) and 
the one with zero proper motion (orbit 2, Fig. [6]) produce 
similar results, whereas the orbit with the minimal orbital 
velocity (orbit 3, Fig. produces a cluster which appears 
largely extended. 



4.1 Absolute magnitude 

From these snapshots we measure the integrated absolute 
magnitude, Mv, of each cluster represent ation withi n a ra - 
dius of 26 arcmin around its centre, just as ICote et al l (|2002T ) 
have done for their observational data. The results are listed 
in Tab. IU3l for the three orbital types. We see that the clus- 
ters starting off with smaller initial masses, independent 
of the orbit, have lost too much mass within the 3 Gyr 
of evolution, such that today their absolute magnitude is 
too low compared to the observational value of Mv = —3.8 
mag. Clusters with Mq 5000 Mq lose just about the right 
amount of mass within this time. From the tables we can ex- 
pect clusters with initial masses even higher than 10000 Mq 
to exceed the observed absolute magnitude. This suggests 
that our range of initial parameters covers the right part of 
the parameter space of initial conditions. 



Table 1. Results for the TV-body comp utations with the orbit 
using the transverse velocity measured bv lSiegel et al.l l l200ll h i.e. 
orbit 1. Mo gives the initial mass of the cluster at the beginning 
of the computations, and Ro its initial half-mass radius. My is 
the measured absolute magnit ude within the in ner 26 arcmin at 
an age of 13 Gyr, i.e. today. ICote et al.l j2002l ) find a value of 
Mv = —3.8 mag for Pal 13. ay gives the velocity dispersion 
within the inner 2 arcmin measured from a sample of 21 stars. 
The value gives the mean of 10 8 independent measurements (see 
Sec. l4.2l for details), the uncertainties give the limits in which 67% 
of all measurements lie. I Cote et al.l j2002t) find a r = 2.2 ± 0.4 for 
Pal 13. rj is the slope of the surface density profile at large radii 
measured with the KKBH template. The uncertainties give the 
standard error from a least square fit. For Pal 13 we measure a 
slope of about 1.9 based on the observational data by Cote et al. 
(2002). 



Mo 
[M ] 



-Ro 

[PC] 



M v (R < 26') o> 



mag 



(R < 2') 
[km/s] 



3000 
3000 
3000 
4000 
4000 
4000 
5000 
5000 
5000 
7500 
7500 
7500 
10000 
10000 
10000 



4.0 
6.0 
8.0 
4.0 
6.0 
8.0 
4.0 
6.0 
8.0 
4.0 
6.0 
8.0 
4.0 
6.0 
8.0 



-2.7 
-2.6 
-2.1 
-3.0 
-3.2 
-2.6 
-4.0 
-3.1 
-4.0 
-3.9 
-4.2 
-3.5 
-3.7 
-4.2 
-3.6 



n =0+0.09 
"■ o - 3 -0.09 
,,+0.09 
-0.09 
-,+0.07 
-0.07 
,+0.11 
-0.11 

63+ 11 
U-tM-O.lO 
ro+0.08 

3 +0.12 
-0.11 
,+0.12 



0.54Z 
0.46j 
0.63" 1 



0.69 4 



0.92^ 



0-711o.il 

U-OS_o.09 
,+0.16 
-0.16 

0-841 1% 
0.73^1? 

1+0.20 
-0.19 
-,+0.16 
-0.15 



1.11. 

0.96j 
0.83 



+0.13 
-0.13 



3.86±0.58 
3.55±0.34 
4.34±0.76 
3.33±0.20 
3.77±0.12 
4.15±0.16 
3.58±0.14 
3.82±0.15 
3.86±0.14 
3.36±0.07 
4.00±0.16 
4.08±0.07 
3.86±0.24 
4.40±0.21 
4.24±0.15 



Moreover, comparing the same clusters but on the dif- 
ferent orbits, we find that the clusters evolve quite simi- 
larly internally and that their absolute magnitudes are only 
marginally influenced by the orbital type. In fact, at the be- 
ginning of the computations the clusters of a given mass and 
size are exactly the same clusters just on different orbits. In 
this way we make sure that differences come from dynami- 
cal evolution and not from stellar evolution. After 3 Gyr the 
masses between the clusters of a given initial mass and size 
differ by only 50 — 100 Mq . From this we can deduce that 
the influence of the pericentre passages on all three orbits 
are rather unimportant. Otherwise, the more eccentric orbits 
(orbit 3 and orbit 1) would have induced more dissolution 
on those clusters, and altered the final absolute magnitudes 
more significantly. 



4.2 Velocity dispersion 

From the computations we also take rad ial velocity disper - 
sion measurements in the same way as ICote et all (|2002l ) 
have done. That is, we draw 21 stars from the sample of 
stars within the inner 2 arcmin of the clusters, while mak- 
ing sure that all 21 stars lie within 10 km/s of each other. A 
star with a radial velocity differing more than 10 km/s from 
the other stars which have been drawn from the population 
would therefore be regarded as a non-cluster member, even 
though this does not necessarily hold true in our computa- 
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Figu re 5. Left: s tellar map of a 4 deg 2 field around Pal 13, computed using the orbit with the transverse velocity measured by 
ISiegel et al.l d200lll . i.e. orbit 1. Each dot represents a star above my = 25 mag. A background of 0.28 stars/arcmin 2 was added with 
random positions. Right: underlying surface density map of Pal 13 for the same field. One bin corresponds to 9 arcmin 2 . The colour 
coding shows log 1 g(A r + 1), where TV is the number of Pal 13 stars in a bin. The expected background of 0.28 [0.69] stars/arcmin 2 
corresponds to a value of 0.5 [0.9] in this representation. At the distance of Pal 13, 1 deg corresponds to roughly 420 pc. The cluster is 
well limited and the density falls off steeply, only small traces of tidal tails can be seen. 
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Figure 6. Left: stellar map of a 4 deg 2 field around Pal 13, computed using the orbit with zero transverse velocity (orbit 2). Each 
dot represents a star above my = 25 mag. A background of 0.28 stars/arcmin 2 was added with random positions. Right: underlying 
surface density map of Pal 13 for the same field. One bin corresponds to 9 arcmin 2 . The colour coding shows log 10 (7V + 1), where N 
is the number of Pal 13 stars in a bin. The expected background of 0.28 [0.69] stars/arcmin 2 corresponds to a value of 0.5 [0.9] in this 
representation. At the distance of Pal 13, 1 deg corresponds to roughly 420 pc. The cluster is also well limited, just as in Fig. [5] Only the 
tidal tails are a bit more pronounced since the cluster moves at a lower velocity and hence the stellar density within the tails is higher. 



tions. The velocity dispersio n is then computed i n the same 
fashion as has been done in iKupper fc Kroupal l|20ld ). We 
independently draw 10 6 sets of 21 stars from each cluster 
and compute for each set the dispersion of the stellar veloc- 
ities. From these values we take the mean, which is given in 
Tab. IU3I The uncertainties of these values give the bounds 
in which lie 67% (la) of all measurements above and below 
the mean. 



From Tab[T]we can see that the clusters on the Siegel 
orbit ( orbit 1 ) yield too low velocity dispersions in compar- 
ison to the observational value of oy = 2.2 ± 0.4 km/s. The 
orbit with zero proper motion yields similar results (orbit 2, 
Tab. [2]). In both sets of computations we achieve the highest 
velocity dispersions of 1.1 ± 0.2 km/s in the most massive 
and most compact cluster of 10000 M© and Ro — 4.0 pc. 
This is expected when we assume that the clusters are in 
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Figure 7. Left: stellar map of a 4 deg 2 field around Pal 13, computed using the orbit with the transverse velocity minimizing Pal 13's 
orbital velocity (orbit 3). Each dot represents a star above my = 25 mag. A background of 0.28 stars/arcmin 2 was added with random 
positions. Right: underlying surface density map of Pal 13 for the same field. One bin corresponds to 9 arcmin 2 . The colour coding 
shows log 10 (A r + 1), where N is the number of Pal 13 stars in a bin. The expected background of 0.28 [0.69] stars/arcmin 2 corresponds 
to a value of 0.5 [0.9] in this representation. At the distance of Pal 13, 1 deg corresponds to roughly 420 pc. The cluster is embedded 
in a far-extending cloud of stars, which originates from the compressed tidal tails getting pushed back into the cluster vicinity as the 
cluster-tail system is being decelerated on its way to apogalacticon. 



Table 2. The same as Tab. [T]but for the orbit with zero trans- 
verse velocity (orbit 2). 



M 
[M ] 



Ro 

[PC] 



M v (R < 26') 
[mag] 



, (R < 2') 
[km / s] 



0.56 
0.92 



,+0.09 
^-0.09 
,+0.09 
-0.09 
1+0.07 

i -0.07 
1+0.11 
-0.11 
,+0.11 
3 -0.11 
,+0.08 
L -0.08 
,+0.12 
-0.12 
3+0-11 
-0.11 
+0.09 
-0.09 
+0.16 
0.16 
+0.13 
-0.13 
+0.11 
-0.11 
+0.19 
-0.19 
+0.15 
-0.15 
+0.12 
-0.13 



virial equilibrium. Note that the observed velocity disper- 
sion is more than 5a off. This fact led Cote et al. to the 
assumption that Pal 13 may contain dark matter or got 
catastrophically heated by the last pericentre passage. 

But for the third kind of orbit, when the cluster is on 
an orbit with a lower orbital velocity such that it is nowa- 
days closer to apogalacticon (orbit 3), the measured veloc- 
ity dispersion is significantly higher (Tab. [3}. We get val- 
ues of about 2.2 km/s within la for many clusters in the 
set. This is due to the number of unbound stars within 
the cluster, so-called potential escapers, and stars outside 
the tidal radius lying in projection within the inner 2 ar- 
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05±0 
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4 


20±0 
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29 
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05±0 


22 
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47±0 


22 



Table 3. The same as Tab.[T]but for the orbit with the transverse 
velocity minimizing Pal 13's orbital velocity (orbit 3). Velocity 
dispersion values which agree within \o with the observed veloc- 
ity dispersion of 2.2 ± 0.4 are bold faced. 



M 
[M ] 



Ro 

[PC] 



My- (R < 26') 
[mag] 



oy (R < 2') 
[km/s] 



'/ 



3000 
3000 
3000 
4000 
4000 
4000 
5000 
5000 
5000 
7500 
7500 
7500 
10000 
10000 
10000 



4.0 
6.0 
8.0 
4.0 
6.0 
8.0 
4.0 
6.0 
S.O 
4.0 
6.0 
8.0 
4.0 
6.0 
S.O 



-2.7 
-2.6 
-2.1 
-3.0 
-3.2 
-2.6 
-4.1 
-3.1 
-4.0 
-3.9 
-4.2 
-3.5 
-3.7 
-4.2 
-3.6 



1.13 
0.48 
+ 0.96 
0.75 
+0.93 
-1.24 
1.06 
0.211 
1.14 
0.47 
+ 1.10 
0.76 
+0.17 
0.17 
+0.39 
0.20 
-1.14 
-0.63 
04+O.I6 
U - M ^-0.17 

90+ 18 
U - MU -0.19 

n sq+1-06 

U - 89 -0.28 

1 16+ 020 
1 - 1D -0.20 

,+0.16 
-0.17 
7+0.56 
-0.19 



0.92 
1.21 
1. 
0.73 
0.99 
1.15 
0.75 
0.80; 
1.11 



0.96Z 
0.87 4 



1.90±0.07 
1.71±0.14 
1.48±0.07 
2.13±0.09 
1.68±0.10 
1.51±0.09 
1.90±0.10 
1.51±0.10 
1.62±0.10 
2.02±0.15 
2.16±0.20 
1.59±0.21 
2.52±0.11 
2.66±0.30 
1.96±0.31 



cmin of the cluster, which pollute the velocity dispersion 
l|Kiipper et al.ll2010bl ). This effect is more significant for the 
clusters which are initially more extended as they show more 
potential escapers. This is a consequence of them being en- 
ergetically more affected by th e pericentre passages (see e.g. 
iGnedin. Lee fc Ostrikenll999l ). 

The clusters on orbit 3 indeed show much more extra- 
tidal material. Looking at the right panels of Fig. I5l7l we see 
that there is barely any stellar material outside the cluster 
for orbit 1 (Fig. [5}, and only little more for orbit 2 (Fig. [6j) . 
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In contrast to this, orbit 3 shows a cluster with an unusual 
extra-tidal extent of several hundred parsec (Fig. [7|). This 
extra-tidal material results from the compression of the tidal 
tails as the cluster and its tails are being decelerated. This 
deceleration is so strong that the whole system consisting of 
cluster, leading tail, and trailing tail, is compressed to a few 
hundred pc. In fact, the system extends even further than 
can be seen in the figure and would extend much further 
if the JV-body computations would have been made for the 
full 13 Gyr since the tidal tails need many Gyr to grow to 
such extent. Note also, that the shape of this system is quite 
irregular since we look at a folded stellar stream and not at 
a bound structure in equilibrium. 

Note that while our choice of the initial density profile 
may well affect the measured velocity dispersion, e.g. in case 
of a higher concentrated King model, it will not affect the 
appearance of the tidal debris. That is, the choice of profile 
may influence the internal structure of the cluster and also 
its mass loss rate but not its tidal debris since the debris is 
formed by orbital com pression and not by the mass loss rate 
(|Kiipper et alj|2010bft . 

4.3 Surface density profiles 

This effect of compression of the tidal debris can also be seen 
as an increase of density within the surface density profiles 
of the modelled clusters, see Figs. 1 81101 The figures show 
the projected stellar number density of the snapshots for all 
clusters except for the most massive ones, meas u red in rings 
around the cluster centres, just as ICote et al.l (|2002T ) have 
done it with their CFHT data (compare with Fig. [!}. For 
better comparison with the CFHT data, we subtract the 
same background of 0.69 stars/arcmin 2 from our data in- 
stead of our lower estimate of 0.28 stars/arcmin 2 since this 
background estimate is important for the outermost data 
points and therefore may influence the fit of the KKBH 
template. Error bars in the figures give the square-root of 
these resulting values as statistical uncertainties. Differences 
in the numbers of stars at small radii within this diagram 
mainly originate from the fact that Cote et al. go down to 
mv = 23.5 mag with their CFHT data whereas we cut at 
mv ~ 25 mag (assuming that Pall3 is observed with an 8m 
class telescope). We fit a KKBH template to these surface 
density profiles in order to measure the slope of the extra- 
tidal material. The results of these fits are displayed in the 
figures (rj, i.e., eta) as well as in Tab. IH3l 

The Siegel orbit (orbit 1) yields well limited clusters 
with steep slopes outside the tidal radius between r\ — 3.3 
and r\ = 4.4 ( Fig. [51), just as expec ted from clusters near 
perigalacticon l|Kiipper et alj|2010rj ). Furthermore, there is 
no clear trend in the slopes with respect to the initial con- 
ditions. The differences are just the statistical fluctuations 
which Kiipper et al. also found in their iV-body data. The 
same holds true for orbit 2 (Fig. 0). Just the scatter is a bit 
larger from r\ = 3.1 to rj = 4.5 but, again, without any clear 
trend. If the slope at large radii was a consequence of the 
last pericentre passage, then we would expect a correlation 
of this slope with the initial half-mass radius, Ro, of the clus- 
ter, since a more extended cluster should be more affected 
by tidal shocking and therefore produce a more pronounced 
tidal debris. 

In contrast to that, orbit 3 results in quite different 
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Figure 12. Comparison of proper motion measuremen t from our 
computations and from the work of lSieeel et al. The data 

points of orbit 1-3 show the proper motions of all stars within 
the inner 6 arcmin of the clusters. As expected, the data points of 
orbit 1 all match precisely with the orbit determined by Siegel et 
al., while the data points of orbit 2 all lie concentrated within the 
origin at zero proper motion. The stars of the cluster on orbit 3 
show a significant intrinsic spread, though, which is due to stars 
belonging to the cluster's tidal debris and which therefore do not 
tightly follow the bulk cluster motion. The data from Siegel et 
al. is shown with and without colour correction to illustrate the 
large spread and uncertainties of the raw data. Also shown is the 
resulting proper motion found by this group. 



surface density profiles (Fig. I10|) . We get a shallow slope at 
large radii for all clusters with values as low as r\ — 1.5, 
and for the steepest not more than 77 = 2.7. That is, the 
clusters near apogalacticon differ significantly from the clus- 
ters which are closer to perigalacticon. Such a behaviour 
of the surface density profile has been found in observa- 
tions of other globular clusters as well. For instance, Palo- 
mar 5, which is kn own to be close to its apogala cticon, 
shows an r\ of 1.5 l|Odenkirchen et al l [200I ). Icii un et al] 
(2010) find 5 Milky Way globular clusters to show shal- 
low values of r\ at large radii, that is 2.44 for NGC 5466, 
1.59 for M 15, 1.58 for M 53, 1.41 for M 30, and 0.62 for 
NGC 50 53. Moreover, AM 4 and Whiting 1 both show a n 
n of 1.8 dCarraro Zinn fc Moni Bid"m1l2007l : ICarrarol[2009l \ 
lOlszewski et al" I (|2009T l furthermore find the Galactic glob- 
ular cluster NGC 1851 to be surrounded by a 500 pc halo 
of stars. Its surface density profile shows a slope of rj = 
1.24 ± 0.66. Our investigation suggests that those clusters 
are all affected by this orbital effect. 



F rom Fig. \T0\ we furthermore find that fitting a King 
l|l962h template to these clusters with the shallow slopes 
at large radii sometimes yields very different results for the 
tidal radius in comparison to the KKBH template (e.g. for 
M = 3000 M Q and R = 8 pc we get Rf KBH = 1.5 ar- 
cmin and R^" 9 = 15.3 arcmin). This we have observed as 
well in the original data (Fig. [T]) and thus would be an ex- 

lanation for the large uncertainties in Pal 13's tidal radius 

Siegel et af]|200ll ; ICote et alj|2002l ). 
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Figure 8. Surface density profiles of the iV-body computations with the transverse velocity measured bv lSiegel et al i.e. orbit 1, 

for all clusters between 3000 Mq and 7 500 Mq. Uncertain ties show the square-root of the number of stars in one bin after subtracting 
a background of 0.69 stars/arcmin 2 like I Cote et al.l (12002) 1 have done (compare with Fig.[T]l. The slopes at large radii as measured with 
the KKBH template are given in the panels i.e., eta). Also given in the panels are the tidal radii as fitted b y the King and the K KBH 
template. All slopes at large radii are quite steep, just as would be expected for a cluster near perigalacticon (Kiipper et al. 2010b). 
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Figure 9. Surface density profiles of the Af-body computations with zero transverse velocity (orbit 2 for all clusters between 3000 Mq 
and 7500 Mq. Uncert ainties show the square-root of the number of stars in one bin after subtracting a background of 0.69 stars/arcmin 2 
like lCote et al.l J2002T> have done (compare with Fig. [TJ. The slopes at large radii as measured with the KKBH template are given in the 
panels (rj, i.e., eta). Also given in the panels are the tidal radii as fitted by the King and the KK BH template. All slopes at large radii 
are quite steep, just as would be expected for a cluster near perigalacticon l lKiipper et al.|[2010bh . 
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Figure 10. Surface density profiles of the TV-body computations with the transverse velocity minimizing Pal 13's orbital velocity (orbit 3) 
for all clusters between 3000 Mq and 7 500 Mq. Uncertain ties show the square-root of the number of stars in one bin after subtracting 
a background of 0.69 stars/arcmin 2 like lCote et al.l (12002) 1 have done (compare with Fig.[T]l. The slopes at large radii as measured with 
the KKBH template are given in the panels (rj, i.e., eta). Also given in the panels are the tidal radii as fitted by the King and the KKBH 
template. In contrast to the other two orbits, this orbit yields surface density profiles with shallow slopes at large radii. This comes from 
the deceleration of the cluster-tail system on its way to apogalacticon and the compression it causes. 
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Figure 11. Upper left panel: overview of the computation results showing the r\ (eta) values measured within the s urface density profiles, 
and the measured radial velocity dispersions, ay, of all models in comparison to values observed bv lCote et al. Each model is 

an independent TV-body realisation of Pal 13 evaluated at an age of 13 Gyr (see Sec. [3j. Other panels: expected velocity dispersions 
after correcting the measurements for 20% (upper right), 40% (lower left) and 60% (lower right) binaries. Adding about 20-40% binaries 
improves the match between the results of orbit 3 and the observations, whereas orbit 1 and orbit 2 necessitate >40% binaries in order 
to agree well with the observed velocity dispersion. 



5 DISCUSSION 

Judging by the different measurements we have performed 
on our artificial clusters (upper left panel of Fig. we see 
that the cluster starting off with Mo = 5000 Mq and Rq — 8 
pc on the orbit with the minimal orbital velocity agr e es bes t 
with the observational values found by ICote et al.l l|2002h . 
The match is not perfect, though, and we do not conclude 
from our findings that Pal 13 has to be on this specific orbit 
with these specific values for mass and radius used here. We 
rather demonstrate that most evidence points to Pal 13 be- 
ing in an orbital phase near apogala cticon (in disagree ment 
to the proper motion estimates by ISieeel et al. I l200lh . As- 
suming this solves the existing problems of Pal 13 without 
further need for dark matter, binaries or changes in the law 
of gravity. 

Nevertheless, although our results reproduce the obser- 
vations of Pal 13 without further ado, we have to check how 
the computational results change when adding binaries to 
the clusters. That is, Pal 13 most likely has a significant bi- 
nary population, which we have not modelled in our A^-body 
calculations, but which should i nflate the velocity dispersion 
further. In fact, iBlecha et al.l (|2004l ) measure significantly 



different radial velocities for some of the Cote et al. sources 
with independent FLAMES data. This might well be due 
to binary motion or, alternatively, to a different quality of 
the instruments used. From 7 stars which are available in 
both studies 2-3 show significant variations, hinting at a 
high binary fraction. Thus, Blecha et al. get velocity disper- 
sion values between 0.6 km/s and 0.9 km/s from samples of 
5-8 stars. 

iKiipper fc Kroupal (|2010h show that for the outer halo 
MW globular cluster Palomar 14 the mean velocity disper- 
sion measured from samples of 17 stars increases from about 
0.5 km/s in the case of Pal 14 being without binaries up to 
almost 4 km/s in the case of it having 100% binaries. Thus, 
adding [20, 40, 60]% binaries to the cluster (i.e. out of 100 
systems in the cluster [20, 40, 60] are binaries) increases the 
velocity dispersion by about [1.0, 1.5, 2.0] km/s when mea- 
sured in this way. Since both, Pal 14 and Pal 13, are quite 
comparable in their expected virial-equilibrium velocity dis- 
persion, the effect of their binary content on observations 
should be comparable, too. In Fig.[TT]we therefore addition- 
ally show the results from our investigation corrected for a 
binary fraction of 20%, 40% and 60%. The corrections, da, 
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of 1.0, 1.5 and 2.0 km/s, respectively, were added to the 
measured velocity dispersions quadratically, i.e., 



Tcorr — \fo% + d(J 2 . 



(10) 



As can be seen in the figure, a binary fraction of 20-40% im- 
proves the match of our models on orbit 3 with the observed 
velocity dispersion, which is a reaso nable binary fraction for 
a globular cluster (|Hut et al.l l"l992'). whereas the models on 
orbit 1 and orbit 2 would prefer values of >40% in order to 
agree with the observational value of the velocity dispersion. 
This would have no effect on the discrepant slope of the sur- 
face density profile at large radii of the models on orbit 1 & 
2, though. 

But why do our proper motion values of orbit 3 differ 
so significantly from the Siegel et al. values? First of all, as 
stated above, we do not argue that Pal 13 has to be on this 
specific orbit but is most likely close to its apogalacticon. A 
family of proper motion values around the values we chose 
may reproduce the observations equally well, since we only 
chose those values in order to minimize Pal 13's orbital veloc- 
ity within the MW such that the orbital compression effect 
is maximal. Moreover, the observational values of Siegel et 
al. may also be influenced by the orbital compression effect 
just as the other measurements. Consequently, their uncer- 
tainties may be largely underestimated. 

In Fig . 1121 w e show the original raw data from 
ISiegel et all |200ll ) which was obtained from photographic 
plates separated by a 40 year baseline. After determining 
the zero point of each photographic plate with 140 potential 
cluster stars lying within the inner 6 arcmin of the clus- 
ter, they identified background galaxies which were then 
used to determine the cluster's proper motion. The spread 
and uncertainties of those 16 background galaxies are large. 
Siegel et al. find that this is due to a colour dependence 
of their proper motion, which is why they apply an ad-hoc 
colour correction to their sample and disregard 2 galaxies for 
which no colour information was available. From this cor- 
rected sample they derive the resulting proper motion and 
the relatively small uncertainties, which are also depicted in 

Fig.m 

In the same figure we also show the proper motion as 
artificially measured for all stars within the inner 6 arcmin of 
the three clusters with M = 5000 M and Ro = 8 pc. The 
stars from the cluster on orbit 1 lie exactly on the proper 
motion values of Siegel et al. (as expected since we used these 
values as input). Accordingly, the stars from the model on 
orbit 2 lie concentrated at zero proper motion. Only the 
stars of the cluster on orbit 3 show a significant intrinsic 
scatter of more than 0.5 mas pc" 1 about the central proper 
motion values of the cluster. This is due to extra-tidal stars 
which are on slightly different orbits than the cluster. Even 
though this scatter does not suffice to explain the discrep- 
ancy it may introduce an additional systematic uncertainty. 
Given that Siegel et al. have not accounted for such a possi- 
ble intrisic scatter, and given the nevertheless large scatter 
in their determined stellar proper motions (see their Fig. 5), 
we argue that in the curious case of Pal 13, the statistical 
and systematic uncertainties in the proper motion determi- 
nation may indeed be (several times) larger than the formal 
error derived by Siegel et al. 

A possible alternative would be that the Galactic po- 
tential which was used in our investigation does not prop- 



erly reflect the true potential of the Milky Way. In or- 
der to keep the proper motion determined by ISiegel et al.l 
l|200lft and move Pal 13 closer to its apogalacticon we would 
have to modify the Galactic potential such that it is sig- 
nificantly stronger at the distance of Pal 13, as this would 
mean that the cluster gets decelerated more strongly and 
cannot get to large Galactocentric distances. But the most 
recent measurements of the Milky- Way potential agree well 
with our choice of potential, finding circular velocities of 
abou t 220 km/s at the distance of the Sun (e.g. iGhez et al.l 
2008), and a possible flattening of the potential of about 
0.87 ijKoposov. Rix fc Hogdl2010h . Such a flattening would 
imply the opposite of what would be necessary for bringing 
Pal 13 closer to its apogalacticon, since the Galactic force 
would be less strong at a given distance from the Galactic 
plane compared to the unflattened case. 



6 SUMMARY AND CONCLUSIONS 

We performed a set of 3 x 15 iV-body computations 
of the low-mass Milky-Way globular cluster Palomar 13 
which shows som e peculiarities in observations. First of all, 
ICote et all \2QQ± measured a velocity dispersion of 2.2 ±0.4 
km/s, which yields a very high mass-to-light ratio of about 
40 due to the cluster's low integrated absolute magnitude 
of only Mv = —3.8 mag. Secondly, Pal 13 shows a shal- 
low slope at large radii within its surface density profile of 
77 = 1.9, making a de t ermination of its tidal radius diffi- 
cult (jSiegel et al.l l200ll : ICote et al.ll2002l ). It has been sug- 
gested in both publications that these effects might be due to 
Pal 13's last pericentre passage on its eccentric orbit a bout 
the Galactic centre. In contrast, iKiipper et all l|2010bl ) find 
by means of iV-body computations that pericentre passages 
barely influence the appearance of a cluster's surface density 
profile. Instead, they find a flattening of the surface density 
profile only for clusters on eccentric orbits which are about 
to reach their apogalacticon, resulting from the compression 
of the tidal debris as it gets decelerated in its orbit. 

We therefore use three different orbits to explicitly test 
these two hypotheses for Pal 13 (Fig.[2|. Fir st, we use the ra- 



dial velocity measured bv lCote et alP i 20021) in combinatio n 
with the proper motion measured by Siegel et al.l (|200lh . 
This yields an orbit in which Pal 13 today is close to its peri- 
galacticon [orbit 1). Secondly, we use only the measured ra- 
dial velocity but set the proper motion to zero. This yields a 
similar orbit which is less eccentric and in which Pal 13 is to- 
day closer to apogalacticon but still not close enough for the 
effect o f tidal debris compression described in lKiipper et al.l 
|2010bh , and which causes the surface density profile at large 
radii to become shallow, to take place [orbit 2). Finally, we 
use the orbit with the proper motion which minimizes the 
orbital velocity of Pal 13, since this yields the orbit in which 
Pal 13 is today closest to apogalacticon [orbit 3). 

As it turns out, the model clusters on orbit 3 can readily 
reproduce Pal 13's peculiarities both in terms of surface den- 
sity profile and velocity dispersion (and thus mass-to-light 
ratio), whereas the two other orbits cannot (Fig. Ill[) . While 
the first two orbits yield clusters with regular equilibrium ve- 
locity dispersion, the last orbit yields an enhanced velocity 
dispersions and a much larger spread in velocity dispersion 
values when measured from a subset of 21 stars. This is due 
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to unbound stars within the cluster (potential escapers), and 
extra-tidal stars which get pushed back into the vicinity of 
the cluster when the clus ter-tail system gets d ecelerated on 
its way to apogalacticon (|Kiipper et alj|2010al lbh. 

With this investigation we would like to stress the im- 
portance of the orbital phase of a cluster on its appear- 
ance. Particularly interesting is the orbital phase just before 
reaching apogalacticon, where the cluster and its tails get 
decelerated and thus compressed such that the stellar den- 
sity, especially in the region around the cluster's tidal radius, 
gets enhanced with unbound stars. These stars can alter the 
slope of the surface density profile at large radii, and increase 
the measured velocity dispersion significantly. 

Since any cluster (or satellite in general) on an eccentric 
orbit about a galaxy spends most of its lifetime close to 
apogalacticon, it is likely to be observed in such a phase. 
It is therefore expected that a good fraction of all satellites 
are affected by this effect of orbital compression of their 
tidal debris. Observations not taking this effect into account 
may therefore assume too large tidal radii and/or ascribe a 
pronounced tidal debris to tidal shocking which in reality is 
only due to the deceleration of the satellite-tail system. And 
in some cases it may even lead to drastic overestimates of 
the dynamical mass, as is demonstrated here for Pal 13. 

Moreover, orbital compression of a satellite's tidal de- 
bris can produce stellar systems, not only star clusters but 
also dwarf galaxies, which may appear largely extended, ir- 
regular in shape and dynamically hot, and thus may be 
misinterpreted, for instance, as bound systems embedded 
in dark matter haloes. Whether this effect can explain, for 
instance, the high mass-to-light ratios of dwarf galaxies such 
as Segue 1, which is currently know n as the darkest ultra- 
faint dwarf galaxy l|Geha et al.ll2009h . has to be checked in a 
future investigation focussing on typical dwarf galaxy orbits 
and morphologies. Since dwarf galaxies are in general more 
diluted than globular clusters we expect that the effect is 
more pron ounced in thos e cas es. Available work i n this di- 
rection bv lKroupal (|l997t ) and lKlessen fe Kroupal (|l998l ) in- 
deed suggests sim ilar issues are releva nt for dwarf spheroidal 
galaxies (see also lKroupa et al.ll2010l ). Segue 1, in fact, has 
recently been found in SDSS data to show a prominent tidal 
debris a nd therefore was re-classified as a dissolving star 
cluster (|Niederste-Ostholt et~aTI I2009I V Its enhanced M/L 
ratio was interpreted as contaminated by stars from the 
Sagittarius stream but may well be due to the orbital com- 
pression effect described here. 

Finally, this investigation poses the question how reli- 
able proper motion measurement for halo satellites are (es- 
pecially if they suffer from the orbital compression effect 
described in this work), or if we u nderstand the potential o f 
the Milky Way correctly, i.e. is the lAllen fe Santillanl (|l99ll ) 
potential which we used in this investigation a sufficient ap- 
proximation? If the observed peculiarities in Pal 13 are in- 
deed due to the orbital phase of the cluster, then either the 
proper motion measurement or the Galactic potential will 
be in question. 



ACKNOWLEDGEMENTS 

The authors would like to thank Patrick Cote for making 
his Palomar 13 data available. Moreover, they are grateful 



to Sverre Aarseth for making his NBODY codes accessible 
to the public. AHWK kindly acknowledges the support of 
an ESO Studentship. 



REFERENCES 

Aarseth S. J., 2003, Gravitational N-Body Simulations, 

Cambridge, UK, Cambridge University Press 
Allen C, Santillan C, 1991, RMxAA, 22, 255 
Angus G. W, 2008, MNRAS, 387, 1481 
Bessell M. S., 1990, PASP, 102, 1181 

Blecha A., Meylan G., North P., Royer F., 2004, A&A, 419, 
553 

Carraro G., Zinn R., Moni Bidin C, 2007, A&A, 466, 181 
Carraro G., 2009, AJ, 137, 3809 
Chun S.-H., et al., 2010, AJ, 139, 606 

Cote P., Djorgovski S. G., Meylan G., Castro S., McCarthy 

J. K., 2002, ApJ, 574, 783 
Dabringhausen J., Hilker M., Kroupa P., 2008, MNRAS, 

386, 864 

Dehnen W., Binney J., 1998, MNRAS, 298, 387 
Elson R. A. W., Fall S. M., Freeman K. C, 1987, ApJ, 323, 
54 

Flower P. J., 1996, ApJ, 469, 355 

Geha M., Willman B., Simon J. D., Strigari L. E., Kirby 
E. N., Law D. R., Strader J., 2009, ApJ, 692, 1464 

Ghez A. M., et al., 2008, ApJ, 689, 1044 

Gilmore G., Wilkinson M. I., Wyse R. F. G., Kleyna J. T., 
Koch A., Evans N. W., Grebel E. K., 2007, ApJ, 663, 948 

Gnedin O. Y., Lee H. M., Ostriker J. P., 1999, ApJ, 522, 
935 

Harris W. E., 1996, AJ, 112, 1487 

Hurley J. R., Pols O. R., Tout C. A., 2000, MNRAS, 315, 
543 

Hurley J. R., Tout C. A., Aarseth S. J., Pols O. R., 2001, 

MNRAS, 323, 630 
Hut, P., et al. 1992, PASP, 104, 981 
King I. R., 1962, AJ, 67, 471 
King I. R., 1966, AJ, 71, 64 
Klessen R. S., Kroupa P., 1998, ApJ, 498, 143 
Koposov S. E., Rix H.-W., Hogg D. W., 2010, ApJ, 712, 

260 

Kroupa P., 1997, New A, 2, 139 
Kroupa P., 2001, MNRAS, 322, 231 

Kroupa P., et al., 2010, A&A, in press. [arXiv: 1006. 16471 
Kiipper A. H. W., Kroupa P., Baumgardt H., Heggie D. C, 

2010a, MNRAS, 401, 105 
Kiipper A. H. W., Kroupa P., Baumgardt H., Heggie D. C, 

2010b, MNRAS, 407, 2241 
Kiipper A. H. W., Kroupa P., 2010, ApJ, 716, 776 
McLaughlin D. E., van der Marel R. P., 2005, ApJS, 161, 

304 

Mieske S., et al., 2008, A&A, 487, 921 

Niederste-Ostholt M., Belokurov V., Evans N. W., Gilmore 

G., Wyse R. F. G., Norris J. E., 2009, MNRAS, 398, 1771 
Odenkirchen M., et al, 2003, AJ, 126, 2385 
Olszewski E. W., Saha A., Knezek P., Subramaniam A., de 

Boer T., Seitzer P., 2009, AJ, 138, 1570 
Robin A. C, Reyle C, Derriere S., Picaud S., 2003, A&A, 

409, 523 



16 A.H.W. Kiipper, S. Mieske and P. Kroupa 



Siegel M. H., Majewski S. R., Cudworth K. M., Takamiya 

M., 2001, AJ, 121, 935 
Simon J. D., Geha M., 2007, ApJ, 670, 313 
Spitzer L., 1987, Dynamical evolution of globular clusters, 

Princeton, NJ, Princeton University Press, 1987, 191 p. 

This paper has been typeset from a TffjX/ WFftK file prepared 
by the author. 



